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Abstract. In this paper we investigate the Proca-field in the framework of Loop 
Quantum Gravity. It turns out that the methods developed there can be applied to the 
symplectically embedded Proca-field, giving a rigorous, consistent, non-perturbative 
quantization of the theory. This can be achieved by introducing a scalar field, which has 
completely different properties than the one used in spontaneous symmetry breaking. 
The analysis of the kernel of the Hamiltonian suggests that the mass term in the 
quantum theory has a different role than in the classical theory. 
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1. Introduction 

Though loop quantum gravity (LQG) is one of the most successful non-perturbative, 
covariant theories for quantum general relativity, there is still the question whether its 
classical limit gives back general relativity (h — > 0) or quantum field theory {G — > 0). 
For this reason we investigate a question that has been answered in quantum field 
theory, namely the case of the massive vector field. In quantum field theory mass is 
generated via spontaneous symmetry breaking due to the fact that the mass term in 
the Lagrangian of the massive vector field (the Proca-field) is non-renormalizable. On 
the other hand, LQG is a non-perturbative theory, therefor the question arises how to 
quantize this theory in the framework of LQG and how are the results related to those 
of the quantum field theory . In this paper we will show that the Proca-field can be 
quantized in a non-perturbative, diffeomorphism covariant way We will see that just 
as in quantum field theory, a scalar field enters the formalism, but this scalar field will 
play a completely different role in LQG. 

The paper is organized as follows. In section 2, we will perform the 3+1 decomposition 
of the Proca-field in curved space-time and in section 3 we will quantize the theory 
using the tools developed in LQG. In section 4 we will discuss the results which will be 
summarized in the end of the paper. 

2. Classical theory 

The action of the Proca-field coupled to gravity has the form 



where Co is the Lagrangian density of the gravitational field, g is the metric-tensor, g is 
its determinant and is a U(l) connection with curvature F_^ b . To apply the framework 
of loop quantum gravity to this system, first we have to do a 3+1 decomposition: 
introduce on the space-time manifold M a smooth function t whose gradient is nowhere 
vanishing and a vector field t a with affine parameter t satisfying t a V a t = 1- This gives 
a foliation of space-time, i. e. each t defines a 3-dimensional hypersurfice E t . Let us 
decompose t a into its normal and tangential part 



where n a is the unit normal of the hypersurfice E f , N is the lapse function, N a is the 
shift vector. Define the induced, positive-definite metric on T, t via 



Now we use (JHJ) and express n a with t a and N, exploit the fact that y/^g = N^fq and 
introduce the pull-backs of the quantities , T>\ to E t respectively A^ = q^Al , V a = 




(1) 



t a = Nn a + N c 



(2) 



q a b = g a b + n a n b 



(3) 
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q h a T>l and define A = t a A^, A l Q = t a A 1 ^. The canonical momenta are the following: 

(4) 

(5) 

(6) 

,so we have primary constraints il , ITq, H N) H a . Putting everything together (for details 
see [2] , [6] , [H] ) we obtain for the Hamiltonian of the Proca-field 



n = 


5S 

SA 


= 


, n* = 


5S 


ITv = 


6S 
5N 


= 


, n a = 


8S 


K = 


SS 
5A\ 




f 6S 





H= I (^NH + N a H a + A^d + A G + 



+^(-^ + ^o^-^(^) 2 )) (7) 
H = -?-tr(2[K a , K b \ - F a6 )[B„, E„\ + ^z(E"E b + B'B b ) + 

+^Q ab A a A b (8) 

H a = F 3 ab E b + e abc E b B c (9) 
G = V a E a (10) 
G t =V a El (11) 

where F a b is the curvature of the SU(2) connection A l a , K a is the extrinsic curvature 
and B_ a is the dual of the Maxwell connection A^. The phase space carries a symplectic 
structure where the (nontrivial) Poisson-brackets are: 

{A\{x),E b {y)} = 5 b a 5)5{x-y) (12) 
{A a {x),E\y)} = 5 b a 5{x-y) (13) 

Since n = etc. must hold, we get the following consistency conditions (secondary 
constraints): 

= {U N ,H} = H + ^^m 2 A 2 :=H (14) 

= {n a , H} = n a + Vqm 2 AA a := H a (15) 
= {n , H} = G- ^fqm 2 A := G (16) 
= {UlH} = G l , (17) 

where we used the notation A = Ao ~^ ~ a ■ One can verify that the above constraints 
have a second class constraint algebra. To show this, first let us introduce the following 
linear combination of the primary constraints: 

n := iTv + in (is) 

fl N := U N + N a U a + A U (19) 

n a :=n a + A a n (20) 
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It is easy to show that the above linear combinations have weakly vanishing Poisson- 
brackets with all constraints, so these are first class constraints. The constraints Gi 
(gravitational Gauss constraint) and lri a + AaQ. + A l a Gi (diffeomorphism constraint) 
are also first class (as in the m=0 case, the first generates infinitesimal SU(2) gauge 
transformations, while the latter generates infinitesimal spacial diffeomorphisms) . There 
are only two second class constraints: G and 7i. Their Poisson-bracket is 

{H(x),G(y)} = 



m 



E a {x)N a {x) ^q{y) 



5(x -y)- ^qJx)m 2 (A a (x) - ^^M- )v^5{x - y) : 



N(y) y/q(x) v aJ v ^ ' \~ w N{x) 
:=M(x,y), (21) 

where means that the derivative should be calculated in the y variable. 
To deal with second class systems, one needs to introduce the so called Dirac-brackets 
instead of the Poisson brackets. In the case of field theories, it is done in the following 
way (see jl] for details): first one calculates the matrix Mij(x,y) := {Ci(x),Cj(y)}, 
where Ci (x) are the second class constraints in the theory. After that one calculates the 
inverse of Mij(x,y) in the following sense (since Mij(x,y) is a distribution): 



d 3 zM ik {x,z){M^) kj {z,y) = 5 l3 5{x - y) (22) 
After this the Dirac-bracket is defined as 

{f,9} D := {f,g} - J dtxdMLCiixMM^&yKCjiy)^} (23) 

In our case, Mij(x, y) is a 2 by 2 matrix with components 
M n (x,y) = M 22 (x,y) = 

M 12 {x,y) = -M 2l {y,x) = M{x,y) (24) 
The inverse of this matrix has the same structure: 
(M^) u (x,y) = (M(-% 2 (x,y) = 

(M^) 12 (x,y) = -(M^) 21 (y,x) = M(x,y), (25) 
where M(x, y) satisfies the following differential equation 



m 



E a {x)N a {x) 



Jq(x)M(x,y) - yfi@)m*U*{x) - ^^- )v^M(x,y) = S(x - y) 



N(y) v ^ v ' v iaJ v ^ v ' \— y ' N(x 

(26) 

thus the Dirac-bracket has the form 

{/, 9}d := {/, 9} + J d 3 xd 3 y ({/, H(x)}M(y, x){G(y),g} - {/, G(x)}M(x, y){H{y),g}) 

(27) 

The above construction has two drawbacks: the first is that the above Lagrangian 
is not gauge invariant. The cause of this is the mass term, since if one replaces m=0, 
we will have a first class constraint algebra. This affects only the U(l) gauge, SU(2) 
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symmetries are still valid (since the mass term contains gravitational variables in the 
form of the metric tensor and its determinant). The second problem is that one is 
lead to a system with second-class constraints. The latter problem is more troublesome 
because the canonical quantization becomes quite difficult due to the fact that it is far 
nontrivial to implement the Dirac-brackets in the quantum theory ([3]). 
There is an elegant way of curing both problems ([3], [I], [5]), and that is to introduce an 
auxiliary scalar-field and modify the Lagrangian to have the following form: 



c m = c G + V=gg ac g bd 

= Cq + Cym + C-m (28) 



Note that the above Lagrangian is gauge invariant if the transformation rule for the 
fields under gauge transformation is 

*4 = 9JA (29) 
5(f) = A (30) 

and the original Lagrangian is obtained via the gauge-fixing d^(p = 0. 

To see that the system is first class, we perform the 3+1 decomposition to the modified 

Lagrangian as well. Using the notations already introduced we obtain 



H m = J^(NH + N a H a + Aid + A G) 

H = -^-tr(2[K a , K b ] - F ab )[E a , E b ] + ^(E a E b + B a B b ) + 



(31) 



7r 2 y/qm 2 



2^/qm 2 



+ ^Ua + da<l>) (A + d^) (32) 



U a = F 3 ab E- + e abc E b B c + d a <P)7t (33) 

G = V a E a - 7T (34) 
Gi = V a El (35) 

where 7r is the conjugate momenta for 0: 

SS Jqm?(A - N"A„ + Crf - N'd a 4>) 
" = ^ = N (36) 

We also have the primary constraints II at = n a = II 2 = Ho = which are the same as 
in the previous case. If we calculate Lljv etc., we find that H,H a ,G b ,G are secondary 
constraints in the theory. These constraints are referred to as the Hamiltonian, the 
diffeomorphism (modulo gauge transformations), the gravitational Gauss and Maxwell 
Gauss-constraints. 

The phase space carries a symplectic structure where the (nontrivial) Poisson-brackets 
are: 

{A\{x),E]{y)} = 5 b a 5)5{x-y) (37) 
{A a {x),E b {y)} = 5 b J{x-y) (38) 
Mx),<Ky)} = 5(x-y) (39) 
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Before we turn to the quantization, it is worth to observe some of the properties of the 
above system: 

- It is easy to check that the above system is first class, i.e. the constraint algebra 
is closed. This is due to the fact that the canonical momenta of the scalar field 
appears in the Gauss constraint. 

- Note that the mass only appears in the scalar constraint, which means that gauge- 
and diffeomorphism symmetries are independent of m. 

- The Hamiltonian is a linear combination of constraints, which is not true for the 
case where there is no scalar field, since there the Hamiltonian is quadratic in the 
Lagrange- multipliers. 

- The scalar field and the Yang-Mills field is only coupled to each other in the scalar 
constraint and only through a derivative term and no scalar mass-term required, 
which means that if we will quantize this system the scalar field will have a totally 
different role than the one introduced via symmetry breaking. 

- The term A a + d a (p is gauge-invariant with respect to the gauge transformations 
generated by ( |34|) . so this is different to the case when we couple a scalar field to a 
gauge field via covariant derivatives (actually its more like an affine field). 

- One cannot replace m = in the Hamiltonian formalism to obtain the usual 
Maxwell-field. This is not unfamiliar in loop quantum gravity, since it resembles 
to the case of the Immirzi parameter (by this analogy we do not mean any deeper 
connection, though). There the connection and the electric field can be rescaled as 
A l a —y A l a + (3K l a ,E l a — > This is a canonical transformation since the Poisson- 
brackets are invariant under this transformation. If we substitute the new quantities 
in the Hamiltonian, we find that the Gauss and diffeomorphism constraints are 
unchanged, but one part of the Hamiltonian constraint will have a term proportional 
with (3 2 + 1. (see eg. [10] . [TT] and references therein). Now consider the following 
canonical transformation: 

ix — > mir <p — > — 
m 

E a ^mE a A a ^^ 
m 

This will remove the m parameter from the mass term and furthermore this 
parameter will appear only in the Hamiltonian constraint, the other constraints 
will be independent of m. 



3. Quantization 

In this section we will summarize the tools necessary to quantize the Proca- field in Loop 
Quantum Gravity. The main advantage of the theory is that it gives a covariant, non- 
perturbative method to quantize any physical system with first class constraint algebra. 
Since the symplectically embedded Proca-field is of this type, we can directly apply the 
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results achieved in Loop Quantum Gravity. From these the most important for us will 
be the quantization of general gauge systems ([7]- [12]) and the scalar field ([T6].[T7]). 
In the next sections we will give a brief summary on how one can apply these results to 
the present case, and we refer the reader to the above articles for details and proofs. 

3.1. Quantization of a gauge field 

Let us consider a Yang-Mills gauge field with a compact gauge group G. The Hilbert- 
space can be constructed in the following way: let 7 be an oriented graph in £ with 
ei,...,eE edges and Vi,...,vy vertices. Let h ei be the holonomy of the G- valued 
connection of the field evaluated along the edge. Let us define a cylindrical function 
with respect to a 7 graph in the following way: 

f^(A):=f(h ei ,...,h eE ) (40) 

where / 7 is a complex valued function mapping from G E . Since the gauge-group G is 
compact, there is a natural measure, the Ashtekar-Lewandowski measure, which enables 
one to have a inner product on the set of cylindrical functions. With this the Hilbert- 
space of the Yang-Mills field is defined as the set of all cylindrical functions which are 
square-integrable with respect to the above measure: 

H:=L 2 (A,dp ALtG ) (41) 

In our case, G = SU{2) x £7(1), so 

*Hg,ym ■= L 2 (A S u{2) , dfisu(2) ) ® L 2 (Au(i) > ^u(i) ) (42) 

In order to analyze the action of the Hamiltonian and to compute its kernel, it is 
convenient to introduce a complete orthonormal basis on the Hilbert-space ( T421 . 
On the space of L 2 (Asu(2), dpsu(2)) these are called spin network functions and defined 
as follows: let 7 G £ be a graph and denote its edges and vertices respectively 
by (ei,...,ejv) and (vi, . . . , Vy). Associate a coloring to each edge defined by a 
set of irreducible representations (ji, . . . of SU(2) (half-integers) and contractors 
(pi, . . . , pv) to the vertices where pi is an intertwiner which maps from the tensor 
product of representations of the incoming edges at the vertex v\ to the tensor product 
of representations of the outgoing edges. A spin network state is then defined as 

N V 

l T (^) >7i,P := (g)ii(M^)) • <g)Pk (43) 

»=1 k=l 

where • stands for contracting at each vertex Vk the upper indices of the matrices 
corresponding to all the incoming edges and the lower indices of the matrices assigned 
to the outgoing edges with all the indices of pk- 

In the case of L 2 (Au(i), dpu(i)) one must simply replace SU(2) with U(l) in the above 
definition - these are called flux network functions ([IS]). Since U(l) is a commutative 
group, we will have the following definition: for each edge of the graph associate an 
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integer n^. Then the flux network function is defined as 

N 

\nA)> 1 x-=\{(K(A)r (44) 
i=i 

What remains is to define the operators corresponding to the connection and the electric 
field on the Hilbert-space. If we want to implement the Poisson-brackets in the quantum 
theory in a diffeomorphism covariant way, we have to use smeared versions of these fields. 
In the case of gauge fields the natural candidates are the holonomy and the electric flux 
respectively: 

h e (A) = V exp J^A (45) 
E(S) = [ *E, (46) 



where e is a path and S is a surface in S. Then the action of the corresponding operators 
will be defined via the following way: 

k{A)f{A) := h e (A)f (47) 
E(S)f(A):=ih{E(S)J(A)} (48) 

Let us make the action of the electric flux operator a bit more explicit for both the 
gravitational and the Yang-Mills part. The latter case is simple, since the U(l) group 
is commutative. On a flux network function the action of this operator is 

N 

E(S)\F >^ H = ih^THK{S,ei)\F > 7 ,„ (49) 
i=i 

where 

{0, ifeflS , = 0orenS' = e modulo the endpoints 
+ 1, if e lies above S 
— 1, if e lies below S 

In the case of the gravitational part introduce the left and right invariant vector fields 
on SU(2): 

Rhf(9) = J t f(z th 9) and L h f(g) = jj{ge th ) 
Using these quantities the action of the flux operator on a cylindrical function is 

E(S)f,(A) = -J2 <S, e)(6 enSMe) R^ + 6 enSJ{e) L^)f y (A) (50) 

e 

where b(e) and f(e) is the beginning- and endpoints of the edge e and is R on the 
copy of SU(2) labelled by e. 
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3.2. Quantization of the scalar field 

Though the scalar field requires more careful treatment, the results are very similar to 
the Yang-Mills case. The main difference is that in this case one has to use so-called 
point holonomies of the form 

U(v) = exp (i<f>(x v )) (51) 
Then a cylindrical function with respect to a 7 graph has the form 

f,{U):=f{U Vl ,...,U Vv ) (52) 
where v±, . . . , v y are the vertices of the graph. Then the Hilbert space will be 

Hu:=L 2 (U,d^) (53) 

Recently it was shown that a more general description should be used because in 
the original approach the configuration variables are periodic functions, which do not 
suffice to separate the points of configuration space (for details see [XT] and [15]). So 
instead of U(v) we define the configuration variable to be 

U(\,v) = exp{i\(f)(x v )). (54) 

Given a graph 7 with vertices (v\, . . . , vn) and real numbers Aj, . . . , Ajv associated to 
the vertices, let Cyl^ be the set of finite linear combinations of the following functions: 

N 

Then the completion of Cyl = U^Cyl^ will be the C* algebra of configuration 
observables. 

Because the point-holonomy is not smeared, the regulated version of the canonical 
momenta 7r is smeared in three dimensions: 

Pb= [ 7T, (56) 
Jb 

where B is an open ball in E. 

The operator corresponding to (152]) is defined the same way as was done in the case of the 
Yang-Mills field, but the operator version of (T56l) needs careful treatment, because the 
functional derivative of U(v) with respect to is meaningless. The precise calculations 
can be found in |16|, the result is 



P B \D(U) > 7)X := -ih XB (v)X(v)\D(U) > J>X , (57) 

where Xb{v) is zero unless v e B and X(v) is the symmetric sum of left and right 
invariant vector fields at U(v). But since we have a Abelian group, we can write that 
X(v) = Xl{v) = Xr(v), so we obtain 

N 

P B \D(U) > jX .= -ihJ^XB^XjlDiU) > 7)X (58) 
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Now let us apply these results to the case of the Proca-field. Since we have an 
577(2) x i7(l) Yang-Mills field and a real scalar field, the Hilbert-space will be 

H = L 2 {A S u(2), dfj lSU (2)) ® L 2 (Au{i), dfMu(i)) ® ^(W, d//t/(i)) (59) 
with basis (later referred to as generalized spin network functions) 

\S > 7 J,p7a= \ T (A) > l3 , p ®\F{A) >^®\D{U) > 7jX (60) 
Remark: We have seen that in the original (not symplectically embedded) case (171) we 
arrived to a second class constraint algebra, thus we had to introduce Dirac-brackets 
instead of Poisson-brackets. Because of this, the above definitions of the momentum 
operators should be modified by replacing the Poisson-brackets with Dirac-brackets. 
Specifically the momentum operators should be redefined in the following way: 

E(S) D f(A) :=fh{E(S)J(A)} D = 

= ih{E(S)J(A)}-ih J d 3 xd 3 y{E(S),C l (x)}(M^% 1 (x,y){C 1 (y)J(A)}, (61) 

where C{ are the second class constraints G and H. This does not modify the properties 
of the momentum operators, since the Dirac-brackets have the same properties as the 
Poisson-brackets. Only the action on spin network functions changes, which becomes 
more complicated since M(x,y) is not explicit and depend also on and E a . In fact 
only the momentum operator of the Maxwell-field changes, since {G, A l a } = {G, E^} = 0. 
The only question is whether the momentum operator defined above is a well defined 
operator on the Hilbert space, since it is not trivial if its action is a cylindrical function. 
The more detailed analysis of this question can be found in section 6. 



4. The Hamiltonian of the Proca-field 



The total (non-smeared) Hamiltonian of the symplectically embedded Proca-field has 
the form 

H M ) (62) 



H = 


[^N(Hge + Hgl + H e + H b 


+ H P 


'Hge - 


= -e abc tr(F ab {A c ,V}) 

K 




'Hgl = 


= ^e abc tr({A a ,K}{A b ,K}{A, 

K 


»V}) 


He = 






Hb = 






H P = 






2y / gm 2 




Hm = 


■^-q ab (A a + d^)(A b + d b c 





where Hge and Hgl are the Euclidean and Lorentzian part of the gravitational, He 
and Hb are the electric and magnetic part of the Yang-Mills Hamiltonian, Hp is 
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the kinetic term and Hm is the mass term (in the following we consider only the 
symplectically embedded Hamiltonian) . To arrive at a well defined, diffeomorphism 
covariant Hamiltonian operator, one first has to get rid of the y/q quantity from the 
denominator. This is achieved via using the identity (introduced by Thiemann): 

-{Al,V} = 2sgn(dete)el (63) 

K 

The next step is to rewrite the Hamiltonian in terms of holonomies and fluxes. To do 
this, one introduces a regularisation scheme, a A triangulation of E. Each tetrahedra 
in A has a base-point v with three outgoing segments Sj i — 1, 2, 3. Denote the arcs 
connecting the endpoints of and Sj by a^, so we can form the loop := Sj o °sj 1 . 
After an appropriate point-splitting, one approximates the quantities in the Hamiltonian 
by their smeared counterparts, substitutes the corresponding operators and takes the 
limit. This can be done with each of the terms appearing in fl62|) . The details of the 
computations are quite lengthy so here we will only present the regulated Hamiltonians 
and refer the reader to [9] or [13] where one can find the regularisation of all terms 
except the mass term. Since the latter is important in our analysis we will sketch the 
derivation of that particular term. The regulated Hamiltonians are 



= t — ^rle ijk tr(h a -.h s Ah~ s \ V]) (64) 

3ihK ^ E(v) v " hl Sk ' u K J 

AgL = ^ES^^^ (65) 

£t 32 N(V) V"^ JKL MNP ^ 

v v ' v(A)=v(A')=v 

X<H L (A)(^)Qs P (A')(^) X 

x [L JK k< MN ~ E(Fjk)1(F' mn )} (66) 



8 »r/ ^^ X ( V )\ 2 



JJK r LMN 



v(A)=v(A')=v 



l.AW , l,A fe , 1. 



XQUA)(V, ^)Ql A A)( V ^ ^ Q sk(A)(v, -) X 

xQ l SLi A<)(v, l)QZ(A<)(v, \)Q^ m (v, \), (67) 

where 

Q k e (v,r) = tr(r k h e [K\V(vY}) 

E(v) = n ("~ 1 K"~ 2 ) ; w here n is the valance of the vertex v 
Fjk is a surface parallel to the face determined by sj and sk 

What remains is the quantization of the mass term. It is easy to see that this term can 
be obtained from the derivative term of (9] with the replacement d<p — > dip + A, so the 
calculation is completely similar. From q y/q = ^ and Ef = e e^-^ we have 

f f 

H M = — I d 3 x / d 3 yN(x)xe(x,y)e ijk e ilm e abc e bef x 
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J(d a <P + A a )ele k c )(x) ((d b( p + A b )eTe])(y) = 

= ^- (^) 4 J N(x)e ijk {d<f> + A){x) A {A j {x), V(x, e) 3/4 } A {A k (x), V(x, e) 3/4 
x f Xe(x,y)e mn (d<P + My)^{A m (y),V(y,e)^}A{A n (y),V(y,ef'} 



Now we shell introduce the familiar triangulation (with v = s(0)) and use the following 
(exploiting the fact that we are dealing with Abelian fields): 

17(1, s(5t)) = exp[t((f)(v) + 5ts a d a (f)(v) + o(5t 2 ))} 

h s (0, 6t) = exp^tfts"^ + o(6t 2 ))} (69) 
With this we have that 

17(1, s(5t))h s (0, 8t)U(l, v)' 1 - 1 = tdts a (d a (f)(v) + A^iv)) + o(5t 2 ), (70) 

so we have 

r {d<f> + A){x) A {A j (x), V{x, e) 3/4 } A {A k (x), V{x, e) 3/4 } « 



A 

« Sm (A))^ m(A) C/(l, f (A)) -1 - l]Qi (A) (^, ^)Qt(A)(^ !j) (71) 

Substituting this into ([BSD an d taking the limit e^Owe have 
. _ m 2 /4x6^ AT(^ ^ ijk il 

« v 7 «(A)=v(A')=v 

x [U(l, s n (A))h Sn(A) U(l, v)- 1 - 1] [17(1, s r ( A'))^ r(A .)f/(l, v)" 1 - 1] x 

XQt(A)(^ |)^,(A)(«, ^l(A<)(^ |)^(A)K \) (72) 

This part of the Hamiltonian looks problematic due to open ends of the holonomies, 
but since U(l) gauge invariance is studied with respect to ( 1341) and {G, A a + d a <p} = 
the term Hm is U(l) gauge invariant (open ends of holonomies are compensated by the 
scalar field, since the latter is defined in the vertices). Further more this term is also 
diffeomorphism and SU(2) gauge invariant, thus during quantization we do not come 
up against any problems. 

The total Hamiltonian of the (symplectically embedded) Proca-field is 

H = H GE + H GL + Hp + H M 

As Thiemann noticed earlier for similar systems, the Hamiltonian is well-defined, i.e. 
it doesn't suffer from UV divergences, and this is achieved not via renormalisation or 
spontaneous symmetry breaking but treating the gravitational field dynamical. 

5. Kernel of the Hamiltonian of the Proca-field 

5.1. Complete solution 



Though this Hamiltonian is quite complicated, there are a lot of relevant informations 
that can be extracted from it. First let's look at the action of the different terms in the 
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Hamiltonian on a generalized spin network state. 

The action of the gravitational term Hq changes the graph (as pointed out in |10| ) in 
a way that it adds additional edges (specifically extraordinary edges) to the graph and 
changes the intertwiners, but it doesn't affect the labels which correspond to the matter 
fields. The other parts of the Hamiltonian describe the matter fields. Their structure is 
similar: they all contain matter operators and the Q operator in some way which encode 
the interaction of the fields with gravity. This operator only changes the intertwiners, 
it doesn't change neither the colorings or the graph itself. The derivative operators - 
the electric part of the Yang-Mills and the kinetic term -, don't change the graphs, only 
the coefficients, but the mass term and the magnetic term does. 

It is not an obvious question whether this Hamiltonian possesses a non-trivial kernel, 
but we will show that the construction of generating a solution to the Hamiltonian 
constraint, which was introduced by Thiemann, can be generalized to the present case. 
Let IT > -r-. ^:= IT >, be a spin-color network state. Then < $| is in the kernel of 
the Hamiltonian of the Proca-field if for all |T > s we have 

< <$>\H\T > s = (73) 

The key observation of Thiemann was that the Hamiltonian of gravity acts as it generates 
so called extraordinary edges (see details in pU] or [U]). and with this the kernel 
can be constructed in the following way: Denote the set of labelled graphs (spin-nets) 
So G (70,^0) which contain no extraordinary edges (these are the "sources"). Then 
compute S n+ i by acting H G on the elements of S n and decomposing them into spin- 
network states. The main advantage of the sets S n that 1.) they are disjoint, i.e. 
S n H S m = 5 mn and 2.) finding a general diffeomorphism invariant solution to the 
Hamiltonian-constraint reduces to finding a solution on a finite subspace. 
Since we are only interested in solutions which are diffeomorphism invariant, we use T[ s ] 
instead of T s , where [s] labels the diffeomorphism invariant distribution. In particular, 
let the ansatz for a solution be of the form 

N 

<v& i : =E E c w <T iw 

Since the Euclidean part of the gravitational Hamiltonian maps from SV> to S n+1 , we 
have that the condition 

N 

E c [s] <T\ [s] H GE \T> W] =0 (74) 
i=i [s]es n i 

is non-trivial if and only if [s'] G [5' n '~ 1 ]. Since Hqe = ^ v N(v)Hge{v) and the above 
equation has to hold for all possible N, we have 

N 

E C M < T\[s]Hge{v)\T > m = (75) 
»=i [s]es n i 

for each choice of finite number of vertices v and spin nets. Thus, we arrived at a finite 
system of liner equations with finite number of coefficients. 
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In the case of the Proca-field we have three fields. Since the orthonormal base is of the 
form \T > <g)\F > <g)\D >, we will first look for analogs of the sets in the case of 
the scalar-field and the Yang-Mills field. 

The case of the scalar field is simple: denote the set S^CU)^) of all colored graphs, 
that all vertices are labelled by zero. Now define S^ n+1 ^(U) by acting with U(l,v) (for 
all possible v) on every element of S^(U). From the simple action of U(l, v) it is clear 
that this is equivalent to that the elements of S^(U) are those colored graphs, for which 
the sum of the vertex colorings are n. If we look at the form of Hm we find that it maps 
from S^(U) to S^(U) U S^ n+1 \U) U S^ n+2 \U). 

The case of the Yang-Mills field is a bit more complicated, because both Hb and Hm 
changes the graph. The former adds two (Yang-Mills) loops with color 1, while the 
latter increases the color of two edges by one (non-existent edges can be treated like 
they were edges with coloring zero). This is why in this case the analogs of S^ n > will 
have two indices. Denote by S^ n,m "' (YM) the set of labelled graphs which have n loops 
with color 1 and the sum of the colors on all edges are m. It is easy to see that 
these sets are disjoint and the action of the Hb and Hm operators are the following: 
while H B maps from S^ m \YM) to S( n+2 > m+2 \Y M) , H M maps from S^ n ' m \YM) to 
S^ m+2 \Y M) U S^ m+1 \YM) U S^ m \YM). It follows from the construction that 
5"(o,o) (jf m ) w iH contain labelled graphs with zero colorings on all edges. Also note that 
the set S'" ,m '(M) is empty unless n > m. 

Now, let the ansatz for the solution to the kernel of the Proca-field be of the form 
<*\-=J2 2 2 q.],W,M< T lH® < < (76) 

hj,k,l [s]GSM lf]eS (m rPk\YM) [d]£SW(U) 

Now (1761) is in the kernel of the Hamiltonian of the Proca-field, if for all [s]', [f]', [d]' 

< V\H\T > [s y ®\F > [f y ®\D > w = (77) 
With the same reasoning as before this condition is non-trivial if 

-[s]' G [S 1 ™ 1-1 ] U [S*"^ 2 ] (the union of the two sets is necessary if one takes the action of 
Hql also into account, since this latter operator adds two extraordinary edges) 
-[d]' e [S^\U)] u [5r(«-i)(?/)] u [S^- 2 \U)] 

-if]' e [5K-2ft-2)(yM)] U \S^^- 2 \YM)} U [^^(yjlf)] U [S^' Pk \YM)] 
5.2. Special solutions 

Since the system of equation (1771) is very complicated, it is useful to take some special 
solutions in order to understand the full theory. Let < 0\ym be the "vacuum" flux 
network state of the Yang-Mills sector, which means that it has no Yang-Mills colors 
on either edge (note that this is not actually the familiar vacuum state as was shown 
in pj)], since we are not dealing with Fock-spaces). Similarly, denote the vacuum dust 
network state and vacuum spin network state by < 0\u and < 0|g respectively. Now it 
is easy to check that the state < ^Ic® < 0\ym ® 0\u is a solution of ( 1771) if < is 
the solution of the gravitational part of the Hamiltonian (this is because these "vacuum 
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states" are annihilated by the corresponding derivative operators in H P and He ym , and 
are orthogonal to every state created by the operators in H Y m and H v ). So these special 
states can be interpreted as pure gravity. 

Now let us check states of the form < 0|g® < ^\ym® < 0\u- It is easy to show that this 
is in the kernel of the Hamiltonian for all < ^Iym' In fact the same is true for states 
of the form < 0| G <g> < ^\ym® < ^\u- These states are obviously nonphysical since 
the expectation value of the volume, area and length operators of these states are all 
zero for all volumes, surfaces and curves respectively. (It is worth mentioning that these 
states are in the kernel of all Hamiltonian which have density weight one and composed 
only from the gravitational, Yang-Mills and scalar fields) 

Let us check whether there are solutions of the form < ^Ig® < ^\ym® < 0|t/, where 
< ^Ig® < ^Iym is in the kernel of Hq + Hym- The answer is yes, if < ^Ig® < ^Iym 
contains only flux networks that have U(l) colors only on the loops, not on the edges, 
since in this case these states are annihilated by the operator H M . These states are in the 
subset of the kernel of the Yang-Mills field coupled to gravity. Since currently we do not 
have a semi-classical description of the above system it will be for future investigations 
to check the physical meaning of these states. But if we look at the limit m — > 0, we 
find that these states will be solutions, since < ^\a® < ^\ym® < 0\uX(v)\<j> >= for 
all |0 > and 

< ^\ G ® < ^\ym® < 0\u(H — Hq — H Y m)\4> >= 

=<^\ g (8)<^\ym® <0\uH M \(l>>^ (78) 



6. Gauge fixing 

We used the symplectically embedded Proca-field to avoid implementing the Dirac- 
brackets in the quantum theory. This approach led to a well defined quantum theory as it 
was shown in the previous sections. This was achieved by introducing an auxiliary scalar 
field to the formalism. The theory we gained is equivalent to the original one since if the 
constraints are solved the scalar field disappears automatically. But this equivalence is 
not manifest if we do not solve the constraints, thus we need to introduce gauge fixing. 
As we shall see this will again lead to a system with second class constraints like in 
the original theory (without the scalar field). Below we outline how the quantization of 
such systems could be handled. 

The key observation is, as we pointed this out in section 2, that the original Lagrangian 
can be obtained via the gauge fixing T>^<f) = 0. The strategy will be to implement this 
condition in the Hamiltonian formalism. If we compare the original Hamiltonian with 
the Hamiltonian of the symplectically embedded Proca-field we find that if we substitute 



d a <P = (79) 
- N 



_ 9 A -N a A , \ 

vr = Vgm 2 AT ( 8 °) 
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into (1521) - (1511) . we obtain the constraints (1H|) - (1TB|) . So if we introduce the two extra 
conditions (constraints) 

C a = d a <P = (81) 

C = vr - v/gm — = 0, (82) 

we arrive to the original case. One may ask how come the original theory have second 
class constraints while the symlectically embedded one has only first class constraints. 
The trick is that, as it was pointed out in [19], the conditions (15TI) and (1821) are also 
constraints and if we include these to the constraint algebra, we obtain a system with 
second class constraints. 

(In the case of second class constraints one must use Dirac-brackets, so it is natural to 
ask what did we gain with the symplectic embedding? Actually, the main advantage of 
this method that we were able to quantize the theory without introducing the Dirac- 
brackets. The Dirac-brackets are only needed when one fixes the gauge.) 
Now what remains is to calculate the Dirac-bracket and implement the two conditions 
( IHTl) and ( |82l) in the quantum theory. First we need the Poisson brackets of the 
new constraints with the existing ones. If we define the same linear combinations 
(|T8|) . (|T9l) . (|20|) as for the original case, we find that these are first class constraints. 
Also with the same reasoning as we did there one finds that Gi and 7i a + GiA\ + GA a 
are also first class constraints. Thus we have four second class constraints: 7i, G_, C a , C. 
The elements of the antisymmetric matrix M-p (x, y) are therefor the following: 



Mg> = {H(x),G{y)} = Q (83) 

M[? = {H(x), C a (y)} = -p$—i>j>)6(x - y) (84) 

171 v Q{ x ) 

M [P = {H(x), C(y)} = -m'v^aOQl, + d a <j>)V a ^5(x - y) + 

+m ^— -— 5{x-y) 85 

y/q(x) N(y) 

m£ ] = {G(x), C a (y)} = -V^5{x - y) (86) 

M% } = {G(x), C(y)} = m 2 ^N a (y)V^5(x - y) (87) 

m£ ] = {C a (x), C{y)} = -V^5(x - y) (88) 



With the inverse matrix (M^)^ 1] (x,y) one can define the Dirac-brackets the similar 
way as in the first chapter, the only difference is that now we have six second class 
constraints and the matrix M-p(x,y) is much more complicated. 

After this we quantize the theory in the following way. The Hilbert space and the 
configuration variables are defined as in the original (first class constraint) case, but we 
have to redefine the momentum operators ^(5*) and Pb (see the Remark at the end of 
section 3). Now let us check how our new constraints can be interpreted in the quantum 
theory. Because of the complicated Dirac-bracket, the precise action of the operator 
version of ( 1521) is left for future studies. The constraint ( 15T|) on the other hand is much 
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more simple. If we look at the regularisation of the mass term, specifically at equations 
fl69|) . we see that 

17(1, s(5t)) - U(l, v) = i6ts a d a (j)(v) + o(5t 2 ), (89) 

so the constraint can be implemented as < \P| is in its kernel only if 

< |(17(1, 6(e)) - Z7(l,/(e)))|V >= (90) 

for all \ip > and all e edge (b(e) is the beginning-, f(e) is the endpoint of the edge). Let 
> be a basis element, that is 

\* >= I 5 >7J,P,a= \ T (A) > 7 ,l?®\F(A) > 7X ®\D(U) > 7iX (91) 

It is obvious that the action of the operator in fl90l) will be 

(17(1, 6(e)) - U(lJ(e)m >= \S >^ XXi -|S >^ x ^ (92) 

where Ai and A2 is obtained by increasing the value of A„ in the appropriate vertex by 
one. Since the condition is implemented in all vertices, we have that the coefficient of 
< SI ?=•;*? in < ^1 is the same as the coefficient of < SI ; ? rv if both \D(U) > ? and 
\D(U) > 7 are elements of S^ n \U) for a fixed value of n. In other words if we substitute 
( 1761) into the above constraint we get the following condition on the coefficients 

c [s],[f],[d\i ~ c W,[/],Ma = ( 93 ) 
for all [s], [/], [d]i, [d] 2 , where [d]i, [d} 2 are in the same set S^(U). This condit ion has 
non-trivial solutions, for example the case where only those coefficients of < S| f pf% 
are not zero where A„ is the same for all vertices. 

This was the case when one first implements gauge fixing, then quantize the system. 
One may ask whether it is possible to first quantize the system, then do gauge fixing. 
This is a problematic issue for the following reasons. Consider the operator versions 
of the constraints 7i, G, C. These, when quantized, are smeared with test functions 
N, A ,A respectively and have the form H = N(v)H v etc. In the case of C a one 
simply uses the operator in (|90~1) . Now consider the operator matrix Mij(v,v') - which 
could be interpreted as the operator version of M\p (x, y) - defined with the help of the 
commutators of the constraint operators: M 12 = [H V ,G V ,] etc. The question is whether 
the inverse of this matrix - defined via the condition J2 V " M ik{v, v")(M)^(v", v') = 
o~ijO~vv> - actually exists. If it does, then with the help of this matrix we can define a 
Dirac commutator - analog of the Dirac-bracket - the following way: 

[d l ,d 2 ] D = [6 h d 2 \- 

~\ E([^'^(^)]( il)/ )^^'^)[Q(^),6 2 ] - [d 2 ,C l (v 1 )](M^\v 1 ,v 2 )[C J (v 2 )Am) 

where the Ci are the constraint operators. What now one has to do is to replace the 
commutators with this Dirac commutator. Further more one has to modify the action 
of the momentum operators. To see why, consider first a dust network state 

JV 

\D{U)> iX .= \[{U{\ l ,v l )) (95) 
i=i 
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which can be interpreted if the operator f| i=1 ([7 (Aj, v *)) acted on the vacuum state |0 >. 
If we have first class constraint algebra we have the following identity: 

N 

P\D(U) > 7jX = [P^UiKvMO > (96) 

i=l 

After gauge fixing this identity is the key to define the new momentum operator Pp: 

JV 

P D \D(U) > 7jX := [P^UiX^v^UO > (97) 
t=i 

The new electric flux operator E(S)d can be defined the same way. The only thing 
one has to do is replace the momentum operators in the constraints and use Dirac 
commutators instead of usual commutators (note that the gravitational momentum 
operator does not change since the Poisson-bracket of E l a is zero with all four second 
class constraints). 

The critical part of the above construction is the existence and uniqueness of 
(M)~- (ti , v'). But even if it would exist, there is the other question whether (jMj) is 
really the operator version of the Dirac-bracket? The answers to these questions are the 
requirement that the above construction works. 

If the above operator exists then the anomalies of the constraint algebra are removed. 
First let us focus on the gravitational variables. Since the gravitational gauge and the 
diffeomorphism constraints are first class, and the momentum operator for the canonical 
momenta E l a does not change after gauge fixing, there will be no gravitational anomalies 
in the theory. Other anomalies will not appear since by construction {Ci,Cj}D = for 
all constraints and if ([Mil is the operator version of the Dirac-bracket, we obtain 

[Ci,Cj\D = [Ci,Cj] — 

-\ E([a,d(^i)](M) fc Y(^i^2)[G,^] - [C j ,C k (v 1 )](M)^(v 1 ,v 2 )[C l ,C i ]) = 

V\,V2 

= [CiA\-l(&A)-[CiA)) = (98) 

so there are no anomalies (the above commutator does not impose a new constraint since 
it is identically zero). This is also true for first class constraints that have vanishing 
Poisson-bracket with all constraints, since then the corresponding operators will have 
zero commutator, thus the previous expression is also zero. The only problem is the 
case of first class constraints that have non-zero Poisson-bracket with the constraints 
since then the structure constants will appear in the Dirac-bracket. This is problematic 
in the quantum theory since the structure constants may become operators which could 
cause anomalies. In all cases factor ordering ambiguities occur becuse we have terms that 
contain the product of three non-trivial operators, but these do not cause inconsistencies 
but only change the results of the theory. 
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7. Mass 

Note that in this theory mass is a parameter, in fact a coupling constant which couples 
the scalar field, the Yang-Mills field and gravity. In the classical Hamiltonian analysis 
one can make the following rescaling: 7r — > tt/to, — > m0, — > A^/m, ^ a — > TnE^, 
which is a canonical transformation. But in the quantum regime, this parameter enters 
the Hamiltonian in a non-trivial way. In this sense it is very similar to the Immirzi 
parameter of the pure gravitational case. In the latter case the Hawking-entropy 
provided a tool that helped fix this parameter ([TD],[TT]), so there is a chance that 
with a similar method one might be able to make predictions on the value of m. 
Another way would be to define propagators in loop quantum gravity, since the poles of 
the propagators could be interpreted as mass. But so far the question of time remains 
unsolved in the theory, leaving this idea for future research. None the less there are 
attempts which could provide a solution of the problem of time, see e.g. [20] and 
references therein. But without further input, mass is an undefined parameter of the 
theory which has to be given from experiments. 

8. Summary and outlook 

In this paper we investigated the Proca-field in the framework of loop quantum gravity. 
It turned out that the tools developed in this theory can be applied to the Proca-field 
if one introduces a scalar field and rewrites the Lagrangian as ( 128]) . But this scalar 
field is quite different than the one used in spontaneous symmetry breaking, since there 
is no need for scalar mass term and it restores symmetry rather breaking it. The 
role of this field is actually to make the constraint algebra first class. After rewriting 
the Lagrangian, the 3+1 decomposition and the quantization was straightforward. The 
resulting Hamilton operator is well defined and diffeomorphism covariant. We provided a 
method to calculate the matrix elements of this Hamiltonian and some special solutions, 
which will may be a starting point on the interpretation of the theory. We showed that 
the parameter m is actually a coupling constant and is very similar to the Immirzi 
parameter. We also showed how to introduce a gauge fixing both at the classical and 
quantum level which gives back the original theory. This method could be useful for 
quantizing general systems with second class constraint algebra. 

Three major questions remain. The first is the question of the mass. In our description 
we treated it as a given parameter and did not obtain any condition which could have 
given it fixed values. Research in this direction could give us values for m which could 
prove useful testing Loop Quantum Gravity. 

The second question is closely related to the first and it is the relationship with theories 
with spontaneous symmetry breaking. This would give some more understanding of the 
Higgs mechanism on one hand and gauge fixing on the other. 

The third is that one should give a more precise and more general method for gauge 
fixing in order to use it in other cases. For example the gauge fixing used in the case 
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of the Proca-field was simple in the sense that it did not affect the gravitational part of 
the theory (it was diffeomorphism invariant etc.). But one could imagine gauges that 
break diffeomorphism invariance where the construction given here cannot be applied. 
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